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Abstract 

We study a moving D-brane in a time-dependent background. There is 
particle production both because of non-trivial cosmological evolution, and by 
closed string emission from the brane that gradually decelerates due to a gain 
^ ■ in mass. The particular model under study is a DO-brane in an SL{2, R)/U(l) 

cosmology - the techniques used extend to other backgrounds. 

1 Introduction 

The study of time-dependent phenomena in string theory has gained momen- 
tum, partly under the influence of impressive improvements in the measure- 
ment of cosmological data (e.g. the detailed mapping of the cosmic microwave 
background). The euclidean framework that underlies perturbative string the- 
ory does not lend itself easily to interpretation in time-dependent settings. It 
would seem that any progress we can make in continuing hard euclidean results 
in string theory into the Lorentzian domain, and in interpreting them sensibly, 
is worthwhile reporting. (See e.g. EQPHHIIIIIEIISI) 

In this letter, we choose to study an exact cosmological solution to string 
theory (see e.g. (ZHHUHl)) exhibiting non-trivial time-dependence in the closed 
string sector, and we supplement it with a non-trivial moving brane that solves 
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the equations of motion of classical open string field theory. Exact results 
in the purely closed string background can be obtained from the study of a 
particular coset conformal field theory, SL(2, R)/U(l), while exact results in 
the open string sector can be obtained from analytically continuing boundary 
states of the euclidean model (see ^n]^J^2]). Two of the challenges we face 
is to identify the correct analytic continuation, as well as to interpret correctly 
the time-dependent physics. 

In section we discuss the exact cosmological background, and in part El 
the particle production due to the time-evolution. In section |1] we set up the 
formalism that will allow us to add the effects due to emission of on-shell closed 
strings from an extra source. Then, in section [51 we discuss the particular D- 
brane source we add, and the closed string emission it generates. We wind up 
with conclusions in section |HI 



2 Cosmological Background 

We study the cosmology with metric and dilaton 

,o ,,/! dudv 

ds = \ka\ 

1 — uv 

e 2 * = ~7~ (1) 
1 — uv 

which is known to correspond to an exact two-dimensional conformal field the- 
ory, the gauged Wess-Zumino-Witten model SL(2, R)k/U(l), where we consider 
negative level k. As such, it can be combined with flat space factors, and the 
euclidean two-dimensional black hole coset SL(2, R)/U(l) to create a critical 
(super)string theory background [7j. 

In the following, we will make good use of a set of alternative coordinate 
systems, which can be defined as follows (we will neglect the overall factor in 
the metric on occasion): 

u = X 1 - X° ; v = X 1 + X° 
, -(dX°) 2 + (dX^) 2 

~ 1 - (Jfl)2 + {Z) 

which is useful to study global properties of the spacetime and exhibits the 
conformally flat nature of the cosmology (up to a non-trivial dilaton factor); 

u = e x sinh t ; v = - 
ds 2 = -dt 2 + tanh 2 tdx 2 



cosh 2 t ' 

which is a useful coordinate system for comparison to the euclidean 2-dimensional 
black hole and 



jT2 i -T2j„,2\ 



u = 


e x t 


V = 


-e~ x t 


ds 2 = 


l + t 2K 






ds 2 = 


l + e 2r i 



{-drf+dx 2 ), (4) 
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Figure 1: The Penrose diagram for the two-dimensional cosmology 



which makes manifest the relation of the cosmology to flat space in Milne-like 
coordinates near the origin t = 0. 

The spacetime Penrose diagram is drawn in figure (1). There are four 
regions separated by horizons at u = and v = 0. Regions I and II are 
free of any singularities; they can be thought of as expanding and collapsing 
cosmologies respectively. Region I is defined by uv < 0, X° > 0, while region 
II by uv < 0, X° < 0. The coordinates defined by equation (3) (and also (4)) 
are suitable to parameterize these regions. In terms of the (x, t) coordinates, 
the horizon corresponds to the surface t = 0, and for the cosmological region I 
< t < co. In regions III and IV (uv > 0) (naively) there are timelike naked 
singularities. No observer in region I can ever meet the singularities; they are 
simply outside the future lightcone of such an observer. However, it is clear 
that the singularities can influence the cosmological evolution in region I. In 
this sense region I is similar to a big bang cosmology. The crucial difference is 
that near t = the spacetime is smooth having no curvature singularity and 
initial data can be defined on this surface. In the (x, t) coordinates, the metric 
near t = is given approximately by the flat metric in Milne-like coordinates 
and the string coupling is constant, while as t — > oo the metric approaches the 
flat Minkowski metric and the string coupling constant vanishes exponentially. 
Thus, asymptotically the cosmological region I becomes identical to a timelike 
linear dilaton background. 

We will also make use of the boost symmetry of the cosmology: (u, v) — > 
(e~ x °u,e x °v). It is easy to see that in the (x,t) coordinates this symmetry 
corresponds to translations in the spatial x-direction. Because of the symmetry, 
the x-momentum is conserved, but due to the time-dependence of the metric 
there is no conservation law for the energy. 

In the following, our idea will be to concentrate on the future region of the 
cosmology, namely region I. This region naturally Wick rotates to the euclidean 
black hole (after changing also the sign of the level), which is a well-studied 
euclidean conformal field theory. We thus concentrate only on this part of the 
cosmology 1 . It should be noted that different attitudes can be taken and have 

1 The other regions, namely (parts of the) regions III and IV, analytically continue to compact 
parafermions, SU(2)/U(l). It would be interesting to investigate how to interpret these regions 
and their parafermionic analytic continuation (and in particular the resulting boundary conditions, 
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been taken toward the proper definition of the conformal field theory corre- 
sponding to this and other Lorentzian backgrounds. One consists in defining 
the coset strictly in terms of the (hyperbolic) quotient of the original group 
manifold (see e.g. ^3] [0]), and a second consists in following closely techniques 
of quantum field theory analysis in curved space for a particular region of the 
Lorentzian space under study. We follow the second method since the eu- 
clidean conformal field theory that results from this analysis is under sufficient 
control to allow for exact results. It is also desirable to study further the di- 
rect Lorentzian WZW model (and its quotient), and to further compare results 
obtained through these two complementary methods. 

3 Cosmological Particle Production 

In this section, we wish to study particle creation in the bulk of the cosmology 
for a scalar closed string field (f) with quadratic effective Lagrangian 

C = -V^ge-^i-g^d^d^-micf) 2 ). (5) 

In the region of interest, namely the cosmological region I, the dilaton field 
can be chosen to be arbitrarily small and so as a first step we ignore terms in 
the effective action describing string interactions. The analysis in this section 
is equivalent to a mini-superspace analysis of strings propagating in the cos- 
mology, and the results obtained here should match the exact string theoretic 
results in the limit of infinite negative level. As we shall argue later, some ex- 
act string theory results at finite level can be obtained by suitably analytically 
continuing the corresponding euclidean black hole CFT results. The resulting 
equation of motion is then given by 

(- -^-^ d^e-^y/^g^du) + m 2 ) = 0. (6) 

As usual the current density (corresponding to a complexified field 4>) 

h = e _2 V=0( W - 0* V) (7) 

is conserved and can be used to define an invariant (modified Klein-Gordon) 
inner product 

(tufa) = -i j y^e-^ifadrfa-d^ifodE? (8) 

(in notation borrowed from J3]). Due to conservation of the current, the 
inner product is independent of the choice of the spacelike slice £ we use to 
compute it. The modes of the scalar field can thus be normalized using this 
inner product, and the resulting norm is independent of time. 

We now turn to solving the equation of motion. In terms of the (x,t) 
coordinates, the wave-operator is given by 

□ = 77^-7(9? + (cotht + tanhi)<9 4 - coth 2 tdl), (9) 
\k\a 

spectrum, etc). We thank Costas Kounnas for instructive discussions on these issues. 
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where we reintroduced the level k of the conformal field theory. It is useful to 
define the parameter y = sinh 2 t. In terms of y, the wave-operator becomes 

The equation of motion can then be written as the following eigenvalue equation 

(□ + m 2 )c/> = 0. (11) 

Because f(y) = e~ , y) satisfies a hyper geometric equation with pa- 

rameters 

a = j + 1 - i\p\/2, b = -j- i\p\/2, c = -i\p\ + 1; z = -y, (12) 
we can easily find the solutions and expand the field in terms of them 

JL d ^ y -M/2 F{j + x _ i\ p \/ 2 , -j - i\p\/2, -i\p\ + 1; -y) 
+d*e- ipx y l ^/ 2 F(j + 1 + i\p\/2, -j + i\p\/2, i\p\ + 1; -y). (13) 

Here 

^ = -^-4j(j + l)/(|fc|a') 

p 2 1 + 4s 2 

|fc|a' (1 ' 

can be interpreted as a mass squared. We assume that we are dealing with 
continuous representations of SL(2,R), or in other words, with delta-function 
normalizable wave-functions, for which j = — \ + is where s > 0. 

Having solved the equation of motion, we turn to defining quantum fields 
at early and late times and obtain the corresponding vacua. This requires to 
identify positive and negative frequency modes at early and late times. 

Early modes 

In the early flat region (in Milne- like coordinates) near t = (and t positive), 
the modes defined by the decomposition in equation (fT3|) are approximately 
given by 

<j) p = de ipx t~ iM 

+d*e- ipx &\. (15) 

This leads us to identify the normalized, positive frequency, early modes 
u% = -^=e ipx y-^/ 2 F(j + l-i\p\/2,-j -i\p\/2,-i\p\ + l;-y) 



satisfying (up,up) = 2n5(p — p'). (The index 'e' refers to early modes, while 
the index T will refer to late modes.) We quantize the field in terms of these 
early modes by defining annihilation and creation operators as follows 



dp 
2^ 



oM + a ? t < • ( 16 ) 



The operators satisfy the usual commutation relations 

[«;,«;,t] = 2n6(p-p'). (17) 

The early vacuum |e) is annihilated by all a^, dp\e) = 0. 
A few comments are in order. 

• As we already remarked, the metric near t = is given approximately by 
the flat metric in Milne-like coordinates. The early modes are positive 
and negative frequency with respect to the conformal time rj, defined by 
\t\ = e 71 near t = 0. Thus the early vacuum is approximately given by the 
conformal Milne vacuum (see H2) in the region near t = 0. 

• The corresponding frequencies near t = are entirely determined by the 
momentum independently of the mass of the field. This fact can be ex- 
plained as follows. The modes correspond to particles whose coordinate 
energy (conjugate to the time t) scales like E 2 — ^ = m 2 at early times. 
Thus as t — > 0, for fixed mass and momentum, the energy is dominated 
by the momentum dependent piece, E e ~ \p\/t, which explains the result. 

• Upon the analytic continuation \p\ — > i\n\, and a change in the sign of the 
level, the positive frequency mode continues to the normalizable mode 
of the euclidean black hole (see e.g. JO]). Thus correlators in the early 
vacuum are expected to be related to the corresponding euclidean CFT 
correlators upon a suitable continuation. 



Late modes 

We define the late-vacuum at time t = +oo as follows. To obtain positive 
and negative frequency modes at late times, it is more convenient to write the 
solutions in the form 



/ 



^be^y-^Fij + 1 + i\p\/2,j + 1 - i\p\/2, 2j + 2; -y' 1 ) 



+b*e-'^yiF(-j - i\p\/2, -j + i\p\/2, -2j; -y' 1 ). (18) 

As t — * +oo, the wavefunctions are given by 

= /^ te ^( 4 )i/24« e -t e -i2rt + b * e -t^( 4 )i/2-<* e -t e +ia s t j ( 19 ) 
J 2tt 

which leads to the identification of the normalized, positive frequency, late times 
mode: 

1 



u 



e ip*y-j-i Fij + x + i |p|/ 2> j + i _ i\p\/2,2j + 2; -y" 1 ). (20) 



This choice of modes defines the late vacuum following the steps described 
above. We make the following observations: 

• In the (x,t) coordinates and at late times, the metric asymptotes to the flat 
Minkowski metric but the dilaton is linear in time. The string coupling 
constant is exponentially suppressed in time. The behavior of the late 
times modes as t — > oo matches onto the mode solutions for a timelike 
linear dilaton background 9 . 

• From the behavior of the modes as t — > oo, we see that the parameter 
s defined by equation (|14j) determines the frequency of the modes or the 
energy of the corresponding particles at late times, u> = 2s. 
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Particle creation 

The early and late times modes can be related via a transformation in the 
arguments of the corresponding hypergeometric functions (notations and con- 
ventions are summarized in the appendix): 

u l p = -Le^y-i- 1 F(j + l+i\p\/2,j + l-i\p\/2,2j + 2;-y- 1 ) 

= 1 Jpx lv i\p\/2 F( -i\p\,2j + 2 

V4S [V 4+i- i |p|/2,j + l-i|p|/2 j 
F(j + 1 + i\p\/2, -j + i\p\/2, +i\p\ + 1; -y) 
v -i\p\/2 T i i\p\ + h~i\p\,2j + 2 

{ 1 - i\p\,l + j + i\p\/2,l + j + i^^' 
F{j + 1 - i\p\/2, -j - i\p\/2, -i\p\ + 1; -y)\ (21) 

leading to the identification of the Bogoljubov coefficients 

a/' = 8(p-p')a(p) 

(3/ = 5(p + p')(3(p) (22) 

with 

a(v ) = M T ( i\p\ + l,-i\ P \,23 + 2 

{P) V2 S 4-i|p|,l+j + i |p|/2,l+i + ib|/2 ; 

R(n) - A /H rf -i\pl2j + 2 

^ ~ V2/ ( l +J ^|p|/2, J + l-^|/2 ) - (23) 
The Bogoljubov transformation can be written as follows 

u\p) = a{p)u e {p) + (3(p)u e (-p)* 

u e {p) = a(p)*u\p)- I3(p)u l {-p)* (24) 

leading to the following relation in terms of the annihilation and creation op- 
Gr£itors 

a l p = a{pTal-(3{pY{a%y. (25) 

As usual the unitarity relation \a\ 2 — |/3| 2 = 1 is satisfied. Since the early 
vacuum is not annihilated by aL it can be thought of as an excited state of the 
late vacuum, and therefore an observer at infinity concludes that particles have 
been created. 

The total number of particles (for a given single field) produced from evolv- 
ing the vacuum at early times to late times is given by 

( C |JV'|e) = 27rf(0) J ^\P(p)\ 2 - (26) 

The infinite factor 27r<5(0) can be thought of as a volume factor and arises be- 
cause of translational invariance in the spatial x-direction. Therefore, f dp\[3{p)\ 2 
determines the number of particles produced per unit volume. In our case, we 
obtain 

2 cosh 2 [7r(|p|/2-s)] 
IPI sinh(27rs)sinh(7r|p|)' 1 ' 

Two limits of phase space parameters are interesting: 
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• In the limit of large momentum (and fixed mass), using the on-shell con- 
dition equation (|14|) . the energy is given approximately by u = 2s ~ \p\ 
and so for large momentum |/3(p)| 2 ~ e~ 2,r ' p ' = e~ 2wu) leading to an expo- 
nential suppression in the production of high momentum particle pairs. 

• In the small momentum limit, the energy is given approximately by the 
mass u = 2s ~ m eff = {\k\a'm 2 /± - l) l l 2 and so \[3{p)\ 2 ~ coth ( ^;( f/2) ■ 
This leads to a divergence in the integral arising from the infrared region. 
For a massive field, the divergence is logarithmic. The surprising result 
is that the infrared divergence persists independently of the mass of the 
field. 

Thus the particle spectrum produced is dominated by low momentum quanta. 
We can understand the result better by computing the rate of pair creation for a 
given mode. S-matrix elements between the early times vacuum and late times 
particle states are characterized by the quantity 7 (or its complex conjugate) 
defined by JH] 

s = _£(PI ^ T(2j + l)T 2 (-j-i\p\/2) 

a(p)* T(-2j-l)T 2 (l+j-i\p\/2)- [ ' 

Thus the absolute value of j(uj) characterizes the rate of creation of a given 
particle pair 2 . This is given by 

2 = coshy(\p\/2-s)] 
m n cosh 2 [vr(H/2 + s)] V ; 

For large momentum, u> = 2s ~ \p\ and the rate is suppressed I7I 2 ~ e~ 2nuj . For 
small momentum though, the rate is of order one: for p ~ 0, I7I 2 ~ 1. Thus 
small momentum particles are readily produced. 

It can be easily seen that when we analytically continue \p\ — ► i\n\, the 
amplitude in equation (|2*5)) becomes identical to the mini-superspace result 
for the reflection amplitude TZ(j, n) of a closed string mode (with quantum 
numbers (j,n)) in the euclidean black hole geometry jlUj . The mini-superspace 
result is equivalent to the k — ► 00 limit of the relevant euclidean CFT two- 
point function. Thus at finite level k, we expect the exact answer for the rate 
amplitude to be obtainable from the euclidean CFT reflection amplitude upon 
the same analytic continuation accompanied also by a change in the sign of the 
level. Similar results have been conjectured to hold in the cases of time-like 
bulk and boundary Liouville theories PQ[21- Assuming that this is true, the 
exact amplitude is given by equation (28) multiplied by a pure phase factor 
given by jTHj 

^ r(i + ifrV) 

v b -o-t (30) 



rn - ib 2 



with v h = r(l - 6 2 )/r(l + b 2 ) and b 2 = -l/(|Jb| + 2). This implies that the 
rate (|29|) determined by the modulus of 7 is exact even at finite level k (except 
for an overall rescaling of the metric due to the level shift). Thus the analysis 
of the particle spectrum we obtained can be carried through beyond the mini- 
superspace analysis at finite (shifted) level k. 



2 Conservation of a;-momentum requires that pairs of particles of equal and opposite momentum 
are produced. 
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The results we found above indicate clearly that cosmological particle pro- 
duction leads to significant backreaction with a breakdown of perturbation 
theory. From the two dimensional field theoretic point of view, the problematic 
region is the infrared small momentum region. The divergence we found in 
the produced particle density can be associated to an instability of the early 
cosmological vacuum. For a massive scalar field, the divergence is milder and 
since low momentum particles are produced at a rate of order one in string 
units 3 , the space quickly fills with such non-relativistic particles. In other 
words, the space rapidly fills with dust-like matter, and presumably it eventu- 
ally collapses. From the string theoretic point of view, the situation is more 
complicated and the resulting backreaction more severe. In this case, we need 
to take into account the dependence of the effective mass on the momenta along 
the extra spatial directions and the Hagedorn density of string states at large 
energy (and oscillator level). We discuss this case next. 



Embedding in string theory 

We saw that for a fixed (cosmological) mass of a given closed string field, and at 
large energy, the rate of particle production is either exponentially suppressed 
as a function of the energy I7I 2 ~ e~ 2lTU3 or of order one I7I 2 ~ 1, depending 
on whether the momentum along the x-direction is large or small. When we 
further add appropriate factors to build a critical string theory, e.g. by adding 
a euclidean black hole conformal field theory (at level \k\ +4) and a twenty-two 
dimensional flat space factor, the rate of particle production at large energy will 
be enhanced due to the exponentially large Hagedorn density of closed string 
states: p{uj) = 6^%, i n this example, Th = 75^2 = and so the 

total rate of particle production p(w)|7| 2 at large energy is exponentially diver- 
gent 4 , for low-momentum modes. For high-momentum modes, the exponential 
suppresion of the amplitude is sufficient to render the particle production rate 
finite, for sufficiently large |fc|, i.e. for sufficiently weakly curved cosmologies. 5 
We will generically run into this phenomenon in cosmological models where 
we can identify a regime in parameter space in which the rate of particle creation 
(per field) is of order one (e.g. for low momentum, high masses), while the 
degeneracy of states (as a function of the square root of the oscillator level) 
grows exponentially with an exponent which is (twice) the inverse of the closed 
string Hagedorn temperature (see e.g. ^S] for the factor of two). The type 
of matter produced in this regime consists of heavy, almost non-relativistic 
strings and is similar to the "tachyon matter" produced during the decay of 
unstable branes |2()| |21][22] . Within a time-scale of order the string length, 
enough particles are produced to lead to large backreaction and a breakdown 
of perturbation theory. Understanding the fate of the cosmology then requires 
understanding strings at very high densities. 

3 The time-scale for the breakdown of perturbation theory is set by the string scale. No matter 
how small the initial value of the string coupling is, once enough particles arc produced, they will 
lead to significant backreaction. 

4 We do not expect the details of the spins and polarizations of the closed string fields along the 
extra directions to alter this conclusion. 

5 We would like to thank Costas Kounnas for correcting an erroneous statement in the first version 
of our paper. 
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The breakdown of perturbation theory indicates that non-perturbative ef- 
fects are important in understanding the fate of these cosmological models. 
One expectation is that such non-perturbative effects will cut-off the exponen- 
tially growing density of string states. Given this, particle production in string 
cosmology can be turned into an advantage in constructing realistic models of 
inflationary cosmology since such particle production may result in re-heating 
the universe [23] HH after an early inflationary phase. 

It is useful (and sometimes appropriate) to think of the inverse of the ex- 
ponent in the rate amplitude as a cosmological temperature. In the example 
we studied, the cosmological temperature is momentum dependent, and for 
small momenta, it is clearly far above the Hagedorn temperature since then the 
cosmological temperature is infinite. Another recent example where a similar 
phenomenon was observed is the < k < 2 cosmology of ^H]- In time-like 
bulk Liouville theory, the situation is slightly better in that the cosmological 
temperature is only double the Hagedorn temperature (still leading to large 
backreaction) [2] . 

String theory seems to naturally produce cosmological solutions in which 
the cosmological temperature is larger than the Hagedorn temperature, lead- 
ing to large backreaction 6 . One can ask whether there are ways to lower the 
cosmological temperature, or to heighten the Hagedorn temperature. We note 
that the cosmological temperature depends on the choice of initial state and 
initial surface. The leading (mini-superspace) behavior may also receive higher 
curvature corrections which at strong curvature could significantly alter the 
cosmological production rate 7 . Heightening the Hagedorn temperature can be 
achieved in non-critical strings |17j (although the linear dilaton behavior would 
need to be short-cut by a Liouville-like potential to avoid strong coupling prob- 
lems). It would be interesting to understand better the generic constraints on 
string theory that lead to these UV catastrophes in cosmology. 

Note that our vacuum generates a lot of particles in the IR regime (at small 
x-momentum p). This indicates that we would need to modify the definition 
of the early cosmological state in the infrared, for a proper time evolution. 
Indeed, as we argued the onset of the instability occurs at an early timescale. 
It would be interesting to understand the physical consequences of such an IR 
regularization at early times. 

Finally, we add a word on methodology. Implicitly, we have taken the 
attitude that a second quantized string theory in this cosmological background 
can be defined, such that we can study a collection of second quantized fields 
with standard methods developed for quantum fields in curved space. This 
should be contrasted with computations in the first quantized approach, in 
which it would only seem to make sense to ask questions about the stringy 
S-matrix. It would be interesting to investigate the formulation of string field 
theory on time-dependent backgrounds, and to see how these two attitudes 
should be reconciled. 

After these remarks on aspects of bulk string cosmology, we turn to the 
study of cosmological effects due to non-perturbative time-dependent D-branes. 

6 We thank Dan Israel for discussions on this topic. 

7 However, as we indicated above, by using the connection between the amplitude for pair creation 
and the exact euclidean reflection amplitude, one can show that no such correction occurs for the 
particular cosmology that we studied - except for the effects due to an overall rescaling of the metric. 
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4 Particle production in the presence of a 
source 

In this section, we develop a formalism to take into account simultaneous par- 
ticle creation due to cosmological evolution on the one hand and due to the 
presence of a time dependent source (in a linearized approximation) on the 
other hand. To derive the relevant particle production amplitude, we follow 
e.g. HEl- 

The most general solution to the wave equation in the presence of the source 
can be written as follows 

-2$, 



x,t) = (f) e (x,t) +i dt'dx'^/^ge'^G R (x,t;x',t')p(x',t') (31) 



where e (x, t) solves the homogeneous, sourceless wave equation and G R (x, t; x't') 
is the retarded propagator satisfying 

(n + m")G R (x,t;x',t') = -i ^ ~ . (32) 

The retarded propagator is given by 

G R (x,t;x',t') = 9(t-t)(e\[if>(z,t),<l>(a/,1?)]\e). (33) 

Here |e) refers to the early vacuum. Therefore, it can be formally computed 
in terms of the Wightman functions G(x,t; x' ,t') = (e\(j)(x, t)cp(x' , t')\e) and 
similarly for G(x' , t'; x, t). In terms of the early modes the Wightman function 
is given explicitly by 

G (x,t; X ',t') = J ^f p{x - x yimU^^ ( 34 ) 

where we defined f\ p \{t) = y~*bl/ 2 _F(j + 1 — «'|p|/2, —j — i\p\/2, —i\p\ + 1; — y) 
with y = sinh 2 1. 

Let us suppose that the source starts and stops. Then if we wait until all 
of the source is in the past, the field solution can be written as follows 



dp 
2^ 



1 

a P + —fcr=rPU(p),P) 
\/2 P 



+ h.c, (35) 
where p(j,p) is the Fourier transform of the source defined by 8 



p(j,p) = / dtdxV^e-™p(t,x)e-^(fUt)y (36) 



and evaluated at energies and momenta satisfying the on shell condition, eq. 
(|14l) . Only the on-shell part of the source contributes to the late quantum 
field. In the case when a localized source persists at very late times, the above 
formula applies only in the limit t — > oo, away from the source. 

Using the Bogoljubov coefficients, we can write the early modes in terms 
of the late field modes. Since u e v = a(p)*u l p — (3(—p)(u l _ p )* , we obtain after 
substitution 

4>{x,t) = [^(( a ; + ^p(p))(a(prn l p -P(-p) U l _;) + h.c).(37) 



3 We summarize completeness/orthogonality relations in appendix 1X1 
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We can now identify the coefficient of u l p as the late time annihilation operator 
including the action of the source to obtain 

a 1 / = (a; + ^p(p))a(pr-(a^-^=p(-pr)P(- P r. (38) 



We see that we recuperate standard results for cosmological particle production 
when the source is zero |13j . and when the early and late modes in the cosmol- 
ogy are identical, we recuperate standard production of particles by a classical 
source ^Sj. Finally, we compute the total number of particles produced: 

(e\N l \e) = 1^ (V|/%)| 2 + ^\p(pW( P )+P*(-pM-pT\ 2 ) • (39) 

5 The brane source 

We wish to study DO-branes in the cosmology in the mini-superspace approxi- 
mation. The effective metric for such a brane is given by 

ds 2 

—5- = dudv, (40) 
Ts 

and so it is flat. Because the effective metric is flat, a DO-brane follows a 
straight line trajectory in the (u, v) plane. One example is the brane at fixed 
X 1 = (u+v)/2 = c. All other brane trajectories can be obtained by performing 
a boost on this brane 

e~ x °u + e x °v = 2c. (41) 

From the point of view of the string frame metric, the branes follow accelerated 
(non-geodesic) trajectories. The reason is that they experience a force from 
the changing dilaton background. In terms of the (x, t) coordinates, the 2- 
parameter family of the brane trajectories can be written as follows 

sinh(x — xq) sinht = c. (42) 

Since such probe branes follow straight lines in the (u, v) plane, they seem to 
be oblivious to (naive) space-time singularities 9 . 

The description of such a brane from the point of view of the cosmological 
observer using the (x, t) coordinates is as follows. For c > 0, the brane enters the 
cosmological region crossing the horizon u = at finite null time vq = 2ce~ x ° . 
In terms of the (x,t) coordinates, this event occurs at t = 0, x — ► 00. The 
initial x- velocity of the brane is infinite 10 (and negative). Subsequently the 
brane decelerates until it localizes at x = xq as t — ► 00. A similar description 
can be obtained when c < 0. The fact that the brane decelerates from the point 
of view of this observer has a natural explanation. The dilaton background field 
is time dependent such that it decreases as t goes from zero to infinity. During 
late times, the dilaton field decreases exponentially and the metric is flat - 
in the (x, t) coordinates the background asymptotically approaches the linear 



9 This behavior is directly related to the fact that DO-branes in the two-dimensional black hole 
pass through singularities without much ado |19j . 

10 The proper velocity of the brane is finite as can be seen from the behavior of the metric near 
t = 0. 
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dilaton background as t — > oo. As a result, the brane gets heavier and heavier 
during cosmological evolution. Conservation of momentum then implies that 
it must decelerate. 

Such a brane leads to a delta-function source for closed string fields localized 
on its trajectory. The source is time dependent and part of it will be on-shell. 
So we naturally expect on-shell emission of closed strings 11 . To obtain the 
precise form of the source, we examine the coupling of the brane to the dilaton 
fluctuations 5& in the linearized approximation. This is given by 

J dtdx^ge- 2 *p(x,t)(5$), (43) 

where 

/ \ S(x — x(t)) , , 

and x(t) denotes the probe trajectory. Similar couplings are assumed for other 
closed string fields as well. To get the amplitude relevant for particle produc- 
tion, we need to Fourier transform the source with respect to the early, or late 
time modes and evaluate it on shell, as described in section 0J 

We proceed to analyze the behavior of the Fourier transform of the source 
both in terms of the early and the late times modes. We have been able to 
compute the source in momentum space for the case c = 0. This family of 
branes is special in that the branes remain localized in the (x, t) coordinates 
and the time dependence results solely from the changing mass of the branes. 
Fourier transforming with respect to the early modes, we obtain 

P*(j,P) = 

+oo />oo 

dx / dt cosh t sinh i<5(sinh(x — xq) sinh t)e ipx 
-oo Jo 

(sinhfr*l p lF(j + 1 - i\p\/2, -j - i\p\/2, -i\p\ + 1; - sinh 2 t) = 

-e ipx °r( -ts,is, -i\p\ + 1 

2 y is-i\p\/2 + l/2,-is-i\p\/2 + l/2' , v ; 

while the Fourier transform with respect to the late modes yields 

P*(j,p) = 

/•oo 

dx / dt cosh t sinh t5(sinh(x — xq) sinh t)e ipx 
Jo 

(smh ty 2j ~ 2 F(j + l + i\p\/2,j + 1 - i\p\/2,2j + 2; - sinlr 2 t) = 

1 t is, 1 + 2is, 1/2 + i\p\/2, 1/2 - i\p\/2 

2 { is + i\p\/2 + l/2,is-i\p\/2 + 1/2,1 + is ) ' { ' 

The phase factor e ipx ° in the amplitudes (|45|) and (J46|) can be simply understood 
as a consequence of translational invariance along the x-direction. The localized 
brane breaks this symmetry spontaneously. For a discussion of the Fourier 
transform of the generic source (and the comparison to the exact one-point 



l-OO 



n We note that both the annulus amplitude and the one point function of a closed string vertex 
operator on the disk are not suppressed by the diminishing string coupling since they are of order 

9l 
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function for Dl-branes in the euclidean cigar geometry), we refer to appendix 



We can now turn to the computation of the number of on-shell particles 
produced by the DO-brane. The generic result in equation (|H9|) contains the 
bulk term, which we have discussed in detail in a previous section, and the 
source term. The latter is relatively suppressed by an extra factor proportional 
to the inverse volume of the spatial x-direction. This is because the source is 
localized in this direction (and, in contrast to the bulk particle production, the 
source production does not occur homogeneously in space). For the number of 
particles produced by the source, we obtain 



The energies and momenta satisfy the on-shell condition equation (|14|). The 
computation can either be done by using the Fourier transform with respect 
to the early modes, and then using the formalism of section |1J or equivalently, 
by using directly the Fourier transform with respect to the late modes (which 
automatically takes into account the cosmological evolution encoded in the 
Bogoljubov coefficients). 

The following comments are in order: 

• In the limit of small ^-momentum, the energy is essentially given by the ef- 
fective cosmological mass, and the amplitude squared becomes |p;| 2 /4s ~ 



the brane source in this regime is of order one, while for massless parti- 
cles we get an infrared divergence. The latter can be associated with the 
static long range fields produced by the source. The effective cosmological 
mass depends both on the momenta along the extra directions and the 
oscillator level of the emitted string state. 

• In the limit of large x-momentum and fixed effective cosmological mass, 
the on shell condition is given by u = 2s ~ \p\ . In this limit we obtain an 
extra exponential suppression factor |p;| 2 /4s ~ cj -2 e -?ra; in the amplitude 
of particle production. Notice that the exponent of the suppression factor 
is only half of that corresponding to bulk pair production in the large 
x-momentum regime. 

• When we embed in string theory, we need to take into account the en- 
hancement factor due to the exponentially large density of states at large 
energy. Now in the case of the brane source, the closed string states that 
are produced are left/right symmetric and their density is controlled by 
the open string density [U . 

The considerations above indicate that the spectrum of closed strings emit- 
ted by the brane consists predominantely of small momentum highly massive 
states. If the brane is wrapped along extra flat directions, conservation of 
momentum along these directions forbids emission of strings with non-zero mo- 
mentum along these extra directions. And for localized branes along flat extra 
directions, the amplitudes (for scalars) are modified by phase factors such as 
e iky ^ w ith the dependence of the production amplitude squared on the mo- 
mentum along these directions arising through the effective cosmological mass. 




(47) 
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As we saw, for higher effective mass there is a power law suppression factor. 
Since then there is an exponentially growing density of closed string states as a 
function of the oscillator level, the spectrum consists mostly of slowly moving, 
highly massive particles. Hence the emitted matter remains localized close to 
the brane. Presumably the mass gained by the brane, after backreaction is 
included, results in the production of such massive strings. For massive string 
states, it is important to consider finite k corrections to the mini-superspace 
analysis we described. We comment on how to embed in string theory further 
below. 

As can be seen from the analysis in appendix B (see equation (|54|l). the 
small x-momentum regime of the c 7^ branes is better in that an exponential 
suppression factor appears for high cosmological mass (and hence energy) given 
by e~ nu '( 1 ~ 2q / 7T \ where the parameter q is defined by cosq = |c|. For c = 0, 
q = 7r/2, and this factor is one, while for |c| < 1, 2q/ir < 1, and an exponential 
suppression factor survives. The reason for this phenomenon can be understood 
intuitively as follows. When c ^ 0, these branes enter the cosmological region 
with initial x-momentum, and the resulting radiation of strings with zero x- 
momentum should now be less favored relatively to the emission of non-zero 
momentum strings, always in comparison for the relative strengths in the c = 
case. For \c\ > 1, the parameter q is imaginary and appears through an 
oscillatory factor in the amplitude squared. For these cases the suppression 
factor seems to be the biggest possible. We note that the trajectories of these 
branes, when continued to the other regions of the cosmology, cross the timelike 
naked singularities. 

Embedding in string theory 

Again, we can embed the model in string theory by adding a euclidean black 
hole factor and extra flat space factors to the target space. It is then straight- 
forward to identify a region of phase space in which the energy grows like the 
oscillator level (with all momenta fixed (and small)), where the particle pro- 
duction amplitude is not suppressed significantly, and where the exponential 
degeneracy of states leads to an exponentially divergent number of particles 
produced. As we pointed out above, the exponential growth of the relevant 
states is governed by the open string density. (One should compare this to the 
less catastrophic behavior observed in |2j for the case of unstable branes.) 

When we perform the exact analysis beyond mini-superspace, we observe 
that an extra factor (see appendix El and |10| ) 



appears in the amplitude for one-closed string production, which leads to an 
extra suppression factor e ~ k + 2 in the particle production rate for k negative. 
Sadly, this will generically be insufficient to lead to finite particle production 
rates. There is a narrow window of opportunity, at very small k (close to 2), 
i.e. a highly curved cosmology of the type discussed in [Ej, in which one could 
swamp the Hagedorn density with this suppression factor. 

Notice again the similarity of the resulting particle spectrum to the "tachyon 
matter" produced during the decay of unstable D-branes in flat space (20]. It 
would be interesting to see if we can identify a tachyonic mode in the effective 




(48) 
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field theory on the brane (compare e.g. to ^S]) an d see if we can map the 
problem to a tachyon condensation problem on the brane as in pj. 



6 Conclusions 



We have argued for an analytic continuation of the euclidean SL(2, R)/U(l) 
conformal field theory and its interpretation as the future part of an interesting 
two-dimensional cosmology. We thus manage to incorporate all a' corrections to 
the relevant amplitudes discussed. We have identified and discussed a generic 
bulk backreaction problem due to the Hagedorn degeneracy of string states. 
The introduction of a non-perturbative source, such as a DO-brane, leads to 
additional particle production, which might be kept under control in highly 
curved backgrounds. The formalism we developed for the interplay between 
a cosmology and a time-dependent source should be useful for other stringy 
cosmological models as well. 

We have touched upon many directions of future research in the bulk of our 
paper. Other regions of the Lorentzian cosmology (and D-branes propagating 
in these regions) can be continued to other euclidean conformal field theories, 
which in turn await an interpretation in the Lorentzian realm (e.g. the Dl- 
branes in the minimal models). The patching together of different "euclidean 
conformal field theory regions" after suitable analytic continuation, and an 
understanding of the matching of D-brane trajectories on the boundaries of 
these regions is an interesting problem. An easier first step might be found in 
the analysis of the boundary conditions induced at the timelike singularities by 
particular euclidean conformal field theories. 

More generically, one can investigate the feasibility of formulating string 
field theory in cosmological backgrounds, the analysis of the early vacuum, and 
the relationship between open and closed string vacua induced by non-trivial 
consistency relationships between bulk and boundary euclidean conformal field 
theories. 
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Formulas 

We introduce the following notation for a product and quotient of Gamma- 
functions 



A Formalities 



r 



tl, t2, ■ ■ ■ , t n 

b\ , 62, • • • , b m 



) 



r(ti)r(t 2 )...r(t w ) 
r(b 1 )T(b 2 )...r(b m y 



(49) 
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We use the transformation rule (i.e. rule for analytic continuation) for hyper- 
geometric functions 

(-z)~ a F(a,a + l-c,a+l-b;z~ 1 ) = 

l-c,a + l-b 

r ( i t — — i )F{a, b, c; z) 

1 — b,a + 1 — c 

c, 1 - c, a + 1 - ft ^wrfc-i^i-cpfc + ! _ 6 + ! _ 2 _ z) (50) 
2 — c, c — o, a 

and we make good use of the integral formulas j^H] 

f°° t?( h \ -t-ij w a + t,b + t,c,-t 

/ F(a,b,c;-z)z dz = T{ — ) (51) 

Jo a,b,c + t 

which is valid when Re(t) < 0, Re(a + 1) > 0, i?e(fo + 1) > and c should not 
be a negative integer, and 

POO 

/ x c_1 (a; + z)~ a F(a, b, c; — x)dx = 
Jo 

^.c,a — c + a,b — c + a.^, , , , , . 

r(- F a -c + a,b-c + a,a + b-c + a;l-z) (52) 

it, a + 6 — c + ct 

which holds when Re(c) > 0, Re(a — c + a) > 0, Re{b — c + a) > and 
\arg(z)\ < ir. 

Completeness 

We assume the following completeness/orthogonality relations for the early 
modes (which should be compared to analogous relations for their euclidean 
counterparts [25]): 

P&p) = / ^p(y,x)e-^(/ | ; | (y))* 
/ ^(e^/|;|(y))(e- ip ' K (^ | (y))*) = |iV(j,p)| 2 (2vr) 2 o(j -f)8(p-p>), (53) 

where we recall that y = sinh 2 i. The factor \N(j,p)\ is a normalization factor 
and the measure ^—ge~ 2 ®dt = dy/2. 

B The generic source in momentum space 

In this appendix, we make some observations that lead to a guess for an expres- 
sion for the source in momentum space for a brane with a generic trajectory. We 
will base our guess on the expression for the one-point function of the Dl-brane 
on the cigar geometry, which is the natural euclidean analog of our DO-brane 
in the cosmology, as can be seen by comparing their trajectories. 

Before turning to this problem, we note an additional partial result. We 
can determine the Fourier transform of the source at zero momentum (p = 0) 



17 



for a generic trajectory (parametrized by c) as follows 



/+oo poo 
dx / dt cosh t sinh t5(sinh(x — xq) sinh t — c)e ipx 
-oo A) 

(sdnht)~ i ' p lF(j + 1 - i\p\/2, -j - i\p\/2, -i\p\ + 1; - sinh 2 t) 

1 cos(2j + l)q 



"Ip=o 



2s sinhvrs 



(54) 



where we defined cosq = \c\ for |c| < 1 and q £ [0,tt/2]. 

With this extra trump in our hand, we turn to the full problem. We note 
that the semi-classical one-point function for the Dl-brane in the cigar ^U] 
can be rewritten, after multiplying the one-point function by the inverse of the 
normalization factor for the euclidean eigenfunctions, as 

<v, „> = r ineo S ™^- is + \ n \/ 2 ) F( N + l \ 

Xdosed ' 2ssinh27rs l l/2 + is + \n\/2, 1/2 - is + \n\/2 } 

' e -r(2j+l) + e inn e r(2j+l)^ ^ (55) 

which at r = can be rewritten, for n even, as 

<Xdose d > = e^» + §) 9 i- r( - M+i , - — ).(56) 

2ssmh7rs ± + zs + H/2, i - zs + \n\/2 

We expect the analytic continuation to map |n| to —i\p\ (since this is the pre- 
scription that maps normalizable euclidean modes to early positive frequency 
modes), and indeed, we find a match with our results for the Fourier transform 
at c = (up to an overall constant). 

To address the generic case, it is important to notice that the Dl-brane on 
the cigar intersects most circular cross sections of the cigar in precisely two 
points, while the DO-brane in the cosmology only intersects a spatial slice at 
one point. Naively, analytically continuing the cigar source to the Lorentzian 
theory would lead to a doubling of D-branes. 

To obtain the correct one-point function for a single source, we propose to 
perform the following rewriting of the one-point function 

< Y i ^> = e +iwn/2 e ine ° - re l n l + 1 ) 

<-Xdosed^ 2ssinh27rs 4/2 + is + |ra|/2,l/2-is + |n|/2 j 

(cosh s(7r + 2ri)e inn + cosh s(tt - 2ri)) . (57) 

We wish to identify the two terms with two possible sources in the Lorentzian 
theory, at c and xq = i(9o + ir) and at — c and xo = iOo respectively. The results 
now match onto the results we obtained for generic |c| = cosq at p = 0, taking 
into account that c = sinr' and q = tt/2 — \r'\, for r' = ir £ [ — tt/2, 7r/2], and 
taking into account just the corresponding single source term. 

Although we haven't been able to do the integral exactly, we thus propose 
an expression for the FT of the generic case (depending on the sign of r'): 

-r(T7 : -i"^^ 1 -^)coshvrs(l^r / -).(58) 



4s sinh its cosh its 1/2 + is — i\p\/2, 1/2 — is — i\p\/2 tt' 

We estimate that the generic integral that should lead to this result can pre- 
sumably be performed by delving more deeply into the representation theory of 
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SL(2, R). We observe that the large energy behavior of the one-point function 
crucially depends on the value of the parameter r'. 

This guess for the full semi-classical one-point function can now be extended 
to the exact one-point function for the DO-brane in the cosmology, by repeat- 
ing the above manipulations on the known jlUj exact result for the one-point 
function of the Dl-brane in the euclidean cigar . This leads to the extra factor 
in the exact result which we quoted in the bulk of our letter. 
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